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PREFACE 


This report covers the work completed on the research project ’’Wing- 
Propeller Interference Studies” during the period January 1 through May 31, 
1983. The work was supported by the NASA/Langley Research Center through 
the Cooperative Agreement NCCl-65. The cooperative agreement was monitored 
by Dr. Chen-Huei Liu of the Low Speed Aerodynamics Division (Mail Stop -2:8^) . 
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LINEARIZED POTENTIAL SOLUTION FOR AN 
AIRFOIL IN NONUNIFORM PARALLEL STREAMS 


By 

R. K. Prabhu^ aiui S. N. Tiwari^ 

SUMMARY 

A small perturbation potential flow theory is ai- piled to the problem of 
determining the chordwise pressure distribution, lift and pitching moment of 
a thin airfoil in the middle of five parallel streams. This theory is 

then extended to the case of an undisturbed stream having a given smooth 
velocity profile. Two typical examples are considered and the results 
obtained are compared with available numerical solutions of Euler' s 
equations. The agreement between these two results is not quite satisfac- 
tory. Possible reasons for the differences are indicated. 


^ Graduate Research Assistant, Department of Mechanical Engineering and 
Mechanics, Old Dominion University, Norfolk, Virginia 23508. 


^Eminent Professor, Department of Mechanical Engineering and Mechanics, Old 
Dominion University, Norfolk, Virginia 23508. 
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unknotHi constants 

distance of the airfoil above or below a surface of 
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airfoil chord 

airfoil lift coefficient 

airfoil pitching moment coefficient 

jet momentum coefficient 

distance representing the width of velocity nonuniformity 
unknown constants 
running indices 

circulation of the image vortices 

width of a stream in idiich the velocity is uniform 

slope of the mean camber line 

perturbation velocity components 

undisturbed free strean velocity 

distance from the airfoil leading edge along the chord 
height of the image vortex above the airfoil chord 
shape of the deflected jet sheet 
angle of attack; 
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0-1 0-1 

vortex distribution representing the airfoil 
circulation around the airfoil 
defined by x * (l-cos9) c/2 
defined by x' ■ (l-cos<^) c/2 
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1 . INTRODUCTION 


The study of aerodynamic characteristics of lifting surfaces in non- 
uniform flow is of considerable practical interest. Wing sections behind a 
propeller experiencing a jet-like velocity profile, and tailplane sections 
of a conventional airplane experiencing a wake-like velocity profile are two 
examples of such problems. Solutions of these problems are complex and re- 
quire simplifying assumptions. Even if the viscous and compressibility 
effects are neglected, the presence of vorticity in the approaching stream 
necessitates the solution of Euler's equations. Being nonlinear, Euler's 
equations require numerical treatment. This has been done by several 
workers, see for example Chow et al. (ref. 1) and Whitfield (ref. 2). 

This problem can be simplified considerably by replacing tie approach- 
ing nonuniform stream by an equivalent stream having a stepped velocity pro- 
file with a finite number of discontinuities. Ihis problem can be solved by 
neglecting viscosity and compressibility and making a potential flow analy- 
sis. Karman (ref. 3) gave the basis for a linearized potential flow analy- 
sis for such problems. Glauert (ref. 4) employed this method to solve the 
problem of an airfoil in the presence of a jet. He replaced the airfoil by 
a single vortex and, by computing the increment in axial velocity and 
streamline curvature at the airfoil, determined the increase in lift of the 
airfoil. Ting and Liu (ref. 5) extended Karman' s method further and 
determined the chordwise pressure distribution on a thin airfoil in a 
nonuniform stream. 

In this note, the basis for linearized potential flow analysis for the 
problem of an airfoil in a nonuniform stream is reviewed. The method of 
solution of the integral equation for the unknown vortex distribution of 
ref. 5 is simplified. The analysis is then extended to cover the case of 
five streams with four surfaces of discontinuity. Finally, the problem of 
an airfoil in a stream of smooth velocity profile is also treated by the 
linearized potential flow analysis. 
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2. FUNDAMENTAL BASIS FOR THE ANALYSIS 


The preeent analysis is an extension of Kansan's method of representing 
the flow past an airfoil in the proximity of a surface of velocity discon- 
tinuity. This method is based on a linearized potential flow analysis. 
Consider two parallel streams with velocities Uq and with the x-axis 
being the undisturbed streamline separating the two streams. Let an airfoil 
be located in the lower stream and let (uq, Vq) and (u^, Vj^) be the perturb- 
ation velocity components in the lower and upper streams respectively. If 
the distorted streamline that separates the two streams makes an angle S 
with the undisturbed streamline, then we have (see figure 1) 

^0 


tan 0 


Uo +UQ Ul +ui 


Retaining only first order terms, we obtain 

— -IL 

Uo " Ul 


( 1 ) 


We also require the static pressure to vary continuously across the surface, 
i.e. , 


i\) + 1/2 p [Uo - (Uq + uq) 


2 

“ Vo] 


Po + 1/2 p [Ui - (Ul + 


2 2 
Ul) - ml 


Again retaining only first order terms, we obtain 


uo Uo - Ul Ul (2) 

These are the two basic conditions that must be satisfied across the undis- 
turbed surface of discontinuity, and form the basis for the entire analysis 
that follows. 

2.1 Airfoil Near a Surface of Discontinuity 

Consider two parallel streams with velocities Uq and Ui with the x-axis 
as the undisturbed streamline separating the two streams. Let an airfoil be 
located at the point P (with OP • a) in the lower stream (see figure 1). 
Glauert (ref. 4) replaced the airfoil by a vortex of unknown circulation F. 
He then showed that equations (1) and (2) are satisfied if the flow in the 
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upper stream is that due to a point vortex of strength (F + K) U^/Uq at P, 
and the flow in the lower streaa is that due to the vortex F at P together 
with its defracted image of strength K at the point P* {see figure 2) 
where 


2 2 2 2 
K - F (Uo - Ui )/(Uo + Ui) - Fa. 


He then computed the increase in the axial velocity and Che streanline cur* 
vature induced by the image vortex and determined the change in lift of Che 
airfoil. 

A logical extension of this approach is to replace the airfoil by a 
vortex distribution Y (x) , 0 < x < c, instead of a single vortex F. Each of 
the vortex elements Y(x)dx forms images as described above. The downwash 
at the mean camber line of the airfoil can then be determined in terms of 
Y (x) and its image ay (x) and, by satisfying Che flow tangency boundary 
condition on Che airfoil, the unknown y(x) can be determined. This 
problem is treated as a particular case of a more general problem of an 
airfoil in a jet of finite width which is considered in Che following 
section. 


2.2 Airfoil in a Jet of Finite Width 

A more interesting problmn is Che flow past an airfoil placed in a jet 
of finite width. Let us consider the general case of three parallel streams 
of velocities U.^ , l]^ , and with two undisturbed surfaces of discontinuity 
(AA & BB) as shown in figure 3. Let an airfoil of chord c be placed in 
the middle stream at a distance a below the undisturbed surface AA. In 
this case the conditions (1) and (2) have to be satiefied at both Che sui'- 
faces AA & BB. Ting and Liu (ref. 5) represented the airfoil by a vortex 
distribution of unknown strength y(x), 0 < x < c, and by satisfying equa- 
tions (1) and (2) repeatedly across the surfaces AA &. BB, obtained an infi- 
nite set of image vortex distributions (see figure 3). The downwash at the 
airfoil chord is then given by 

v(x) - I (06) j [ 

0 j.Q (x-x* )^+4( jh+a) ^ (x-x* )^+4( jh+b)^ 


5 


+ 


2 bl B(x-x*) 
(x-x')^ +4( j+1)^ 


h2 


■]} Y(x')dx' 


ORIGINAL fag;i m 
OF POOR QUALITY 


The flow tangency condition at tn^ airfoil requirea 


v(x) ■ IJo (® - m(x)) 


(3) 


(4) 


where a is the angle of attack and m(x) is the slope of the mean camber 
line of the airfoil. Equations (3) and (4) form an integral equation for 
the unknown y(x) subject to the boundary condition (Kutta condition) y(c) 
■ 0. Note in equation (3) the first tern in the integrand is the familiar 
singular term that appears in the classical airfoil theory; the other terms 
are not singular. 

As the first step in the solution of (3), x and x' are replaced by 
6 and ^ using the following transformation: 

X ■ (1-cos 9)c/2; x' ■ (l-cos4)c/2 


Equations (3) and (4) together then result in 


a - m(9 ) 



{ 


cos^-cos6 


(a6)j [ a(cos^cos9) 

(cos^-cos 9)2+16( jh+a)2/c2 


B (cos^-coc8 ) 

(cos^-co80 +16( jh+b)^ /c^ 


+ 


2aB(cos4-cos9) 

(cos4-cos 6)2+16( j+l)h^/c^ 


]} 


y( 4) sin 4> d^ 


(5) 


Ting and Liu (ref. 3) then expressed y(^) as the following series: 


Y(«) 


2 Uo (Ao cot ^/2 * 


m 

I 

m“l 


Am sin in^) 


( 6 ) 
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M 

i 

f 

t 
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whfire A , ffl ■ 0,1,2,...,M are the unknown conatanti. Sinilarly they 
n 

expanded the expresaion a - n(6) into a jrier aeriea aa followa: 

m 

a - «(0 ) “ dn I d coa i» 0 (7) 

• ® 
m"l 

where d , m 0, 1,2, . . . ai.e the known conatanta which can be deterained 

B 

for given a - a(6). Theae expanaiona were then aubatituted in equation 

(^). The right aide of the reaulting equation waa reaolved into a coaine 

Fourier aeriea and, by equating the coefficienta of like Fourier coioponenta, 

a aet of algebraic equationa for the unknowna A a ■ 0,1,2, waa 

obtained. The aolution of theae equationa gave the valuea of A and hence 

in 

the diatribution y('^). 

Thia method of aolution ia rather lengthy. Since the integrand in 
equation (5) ia no more aingular than the integrand in the cla^aical airfoil 
theory, all the methoda of aolvxng the claaaical airfoil integral trould be 
applicable in the preaent caae alao. In particular diacretization of y(^) 
would be poaaible. 

We therefore uae Lan' a method (ref. 6) and diacretize y( 4>). In thia 

method vortex pointa and the control pointa 3. are choaen aa 

k 1 

followa: 


(2k-l) X/2N, k - 1,2...,N 

(8) 

i x/N, i-l,2...,N 

(9) 


With thia, aquation (5) reducea to 
o ~ a(3 , ) - { 


, • a(coa^, -coa8 . ) 

1 r ^ Q \ j r k I 

-coa0 . * Ucoa^ -coa0 .)^+16( j+l)h^/c‘ 

1 j“0 k 1 

2o6(coa^ ,-coa0.) 

k I k ^ . 1 1 

* (coa^j^-coa0 +16( jh+b)^ /c^ * (coa^^^-coaS ^)^+16( j+1 )h^/c^ 


2N ' coa''*., . , - 

k*l k 1 j“0 

B (coa^ , -coa0 .) 
k 


Y(^j^) ain i-1, 2,.. ., M 


( 10 ) 


Thia ia a aet of linear aiaultaneoua equationa for the unknowna Y(t^)» k * 
1,2,. ..,N which can be solved eaaily. The lift and pitching moment of the 
airfoil are then given by 


f 


4 
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dl 

2N 


k-1 

M 

I 

k-1 


c ■ -rr I y(^J *in ♦. (l“Co« ♦.) 


( 11 ) 


( 12 ) 


2.3 Airfoil in Che Middle of Five SCreama 

Following GlauerC (ref. 4) let ua coaaider Che problem of an airfoil in 

Che presence of an infinlce aeries of jeca of che same widch b. Lee che 

velociCy in Che nCh jeC be denoced by U , < n < -. Let an airfoil 

n 

represenCed by a vorCex of circulation T be located on Che axis of the 
principal jet in which Che velocity is Uq . Uten Che flow in any jet can be 
represented by an infinite column of equiapaced point vortices at the 
centers of the jets. Flow in the nth jet is due to the vortex system 
as shown in figure 4. In general 


"k - 0, n 0 (13) 

n ’ 

°Kq - r. (14) 


Now by applying conditions (1) and (2) at the surface of discontinuity 
between the nth and the (n'*’l)th jets, the follcwing fundasental relation for 
the strength of vortices can be obts'>vd. 


/(Itx^ ) 

n+1 


n+1. 


n+s+1 


n+a+1 


- a 


n+1 


n-8 


- < n < - 

- < s < - 


( 15 ) 


This equation can be solved; but as pointed out by Glauert, the complete 
solution is very coaq>lex. If, however, we consider only 5 screams as shown 
in figure S, the problem gets simplified to some extent. Hhen the condi- 
tions (1) and (2) are applied to the four surfaces of discontinuity, we 

obtain a relation similar to relation (IS), except that in this case n 

2 1 0 -1 

assumes the values -2, -1, 0 and 1 >nly. Denoting K , K , K , K and 
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-2 

K by I , J , r , K and L respectively, equation (15) for the 'our dif- 

S S S 8 8 8 

ferent values of n can be written as follows. 


*2 *8.2 ■ ■'s.2 " “2 •’.I-8 


Bo's - •'s - “0 >'-.-1 


B -> Vl - Vl -"-! ^-,-2 


•hereo • (l^ ,)/(0^ +U ; 8 - /(l-a^) n«2, 

n a n+1 n n+1 n n 


I - 0, s > 2. 
s 


Jl - 0, To « r , K_1 =0, 


and L * 0, s < -2. 
s 

The solution of these equations can be obtained by substituting successive** 
ly, positive and negative values of s. However, with some algebraic manip- 
ulations, it is possible to obtain recurrence relations for F ; see 

s 

Appendix A for details. These relations are: 

^8 ’ ^1-s T3_g - ® ^ 


' “0 r.j.g - a_i r_3_g - ao 0-1 ^ "2. 


with To “ r , Ti “ oiT , and T-i * ao F . 


Some of the values of F^ computed using these relations are given below: 


F2 “ - (ooa-i )F = F_2 


F3 ■ (026^ ~ ao<^^ , 
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2 2 2 2 2 
(oQai - aoa 20 i -a|a_iBo^^ “ ^_4, etc. 


la general F * a F where a 's are constants. Now, following the procedure 
n n n 

adopted in the previous section we replace the airfoil by a vortex distribu- 
tion y(x), 0 < X < c rather than by a single vortex. Ihen the defracted 

images will also be vortex distributions a y(x). We note that it is not 

n 

necessary to place the airfoil on the axis of the central stream, the 
images formed when the airfoil is offset from the centerline will be shown 
in figure 5. 

The downwash at the airfoil is then given by 


a (x-x’) 




(19) 


n ^ 0 


Tiransforming x and x' into 9 and respectively as before, we 
obtain 


TT n * * a (cosi|>-cos0) 

^ fcos^-co.!ey^. T 6 ' ii^h^rc ^ ^ YU.)sin* d4. (20) 

^ n a — OD 


nt 0 


On discretizing y(<^), and choosing <t> as the vortex points and 6. as 

iC X 

the control points, as in equations (8) and (9) we obtain: 


n » oo 


j .. j .. a^(cos(^j^- cos 9^) 

^^®i^ 2N ^ ^co8<^ -cosQ . * ^ (co8<|> -cos0 .)^+16n^h'^/c^^^^*^k^®^^^k 

k»l k in k 1 

n ^ (j 


i » 1,2,.. .,N 


( 21 ) 


For the linearized boundary condition on the airfoil at the control points, 
we have 


v(0 p - Ub (<* - m(9 .)), i - 1,2,. . 


( 22 ) 
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Now this set of linear simultaneous equations can be solved for the unknowns 
k ■ 1,2,...,N. The lift and pitching moment coefficients can then be 
computed using the relations (11) and (12) respectively. 

2.4 Airfoil in a Stremn of Smooth Velocity Profile 

So far we have considered an airfoil in a stream with a stepped 
velocity profile. Vfe can extend our analysis to the case of a stream with a 
smooth velocity profile. Let us consider, as in the previous section, an 
infinite series of jets of the smse width h, the velocity in the nth jet 
being denoted by U^. In the limit as h tends to zero, the velocity 
profile tends to the given smooth velocity profile. Let the airfoil be 
placed on the axis of the primary jet in which the velocity is Ug (see 
Figure 6) . The strength of image vortices is then given by the relation 
(15). Glauert attempted to extend this to the case of a stream with 
continuous variation of velocity; unfortunately, since he had only a single 
vortex representing the airfoil, he ended up with a simple logarithmic 
singularity in his expression for the increase in axial velocity as well as 
for the streamline curvature at the airfoil. However, if we represent the 
airfoil by a vortex distribution as was done in the previous sections, and 
satisfy the flow tangency boundary condition at the airfoil, we do not 
encounter any difficulty. This has been done in the following. 

As we noted earlier, the solution cf equation (15) for the case of an 
infinite set of jets of the same width h is very complex. However, if we 
assume small variation of velocity from jet to jet, we can write as 

- if , 2U u + u^ u 

a - -il ELS S ( 23 ) 

n if+lf, 2lf+2Uu+u2 U 

n n+1 n n n n n 

where u ■ (U - U ) << U . 

n n+1 n n 

With Of « 1, we can obtain a first order solution (consistent with the 
n 

linearization done elsewhere in the analysis) for equation (15). The re- 
sulting image system for the primary stremn is then found to be 


11 


Fq - r , at y - 0. 
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F, , y » (2n-l)h, n > 1 

zn~*i n n 


a F , at y ■ (2n-l)h, n < 0. 
n n 


Replacing the airfoil by a vortex distribution y(x) instead of a single 
vortex F , we end up with a similar image system. The downwash at the 
airfoil with this image vortex system is then given by 


v(x) 


a (x-x') 
n 


2n ^0 (x-x* )^+y^ 

n*l,3,... •'n 


a (x-x') 




(24) 


This is the governing equation for the downwash velocity v(x) at a point 
X on the airfoil chord in terms of the unknown vortex distribution y(x) . 
For a given airfoil at a given (small) angle of attack, the slope of the 
mean camber line is known. Then, by satisfying the flows tangency condi- 
tion on the camber line, the unknown y(x), 0 < x < c can be determined. 


The stepped velocity profile is only an approximation to a smooth 

velocity profile. We can formally proceed to the limit of a smooth velocity 

profile by increasing the number of steps (N) and correspondingly decreasing 

the width of the step (h) . For small h (*dy) , u^ ■ (dU/dy) dy, and the 

expression (23) for a becomes 

n 


1 dU ^ 

U dy ^ 


(25) 


where U and dU/dy are measured at (2n-l) h/2 ■ ndy - dy/2, and the corre- 
sponding image is located at (2n-l) h ■ 2ndy - dy. In the limit as N 
tends to infinity, the summations in the integrand in equation (24) are re 
placed by the corresponding integrals. With this, the downwash equation 
(24) can be rewritten as 
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v(x) 


1 f 1 r 1 dU (x~x* )dy 

2it Jq ^ x-x' ” •’o U ^ (x-x')‘^+4y^ 
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+ 



liSL 

U dy 




(26) 


It can be shown that for U and U' of the order of unity, the linearized 
boundary condition is 


v(x)/U(0) •a - m(x) 


(26a) 


For any given smooth velocity profile U(y) with U(y)^0, -«> < y < «, the 
integrals within the brackets can be evaluated using any standard technique 
and the vortex distribution y(x) can be determined by the usual procedure. 
Relations (11) and (12) can then be used to determine the lift and pitching 
moment coefficients. 

If the given velocity profile U(y) is symmetric and the airfoil is 
placed on the line of symmetry, equation (26) reduces to 


v(x) 


J, 

2n 



2(x-x') 


r ^ ^ dy 

^ (x-x' )^+4y^ 


} Y(x' )dx' 


(2 7) 


The unknown vortex distribution can now be determined by the usual procedure 
after which relations (11) and (12) yield the lift and pitching moment co- 
efficients . 

If the given velocity profile U(y) is synmetric and the airfoil is 
placed on the line of symmetry, equation (26) reduces to 


v(x) 


__1 

2n 



2(x-x’ ) 


r iiiL dy 

» “ '*y (x-x')2*4y2 


} Y(x')dx' 


(2 7) 


The unknown vortex distribution can now be determined by the usual procedure 
after which relations (11) and (12) yield the lift and pitching moment co- 
efficients. 


2.5 Jet-Flapped Airfoil in a Nonuniform Stream 
Let us consider a thin jet flapped airfoil in the middle of five 
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parallel sCreama. Following Che n»chod of Spence (ref. 7) the airfoil is 
represented by a distribution of vortex with density y(x) , 0 < x < c, and 
Che jet emerging from Che trailing edge of Che airfoil is represented by a 
distribution of vortex with density Yj(x), c < x < •. As in the previous 
examples, this vortex distribution also forms a system of images. Ihen Che 
expression for the downwash on Che x-axis can be written as follows: 

■ ar ^0 f ^ " jx-x ' )^.A.V 1 

n ■ r “ 

n ^ 0 

= if 0~ 


The coefficients a *8 are Chose described in section 2.3. The boundary 

n 

conditions are: 

on the airfoil v(x) ■ Uq [® - m(x)) , 0 ^ x < c, (29) 

aiul on the jet v(x) ■ -Uq y'(x) , c < x < *. (30) 

Now Yj i-S related Co Che curvature of the jet by the relation 

Y .(x) - -1/2 UoC.y" (31) 

J 

where C is the jet momentian coefficient and y(x) is the shape of the 
%) 

jet. Combining these, we obtain the following pair of integro-dif ferencial 
equations: 

c C * 

/ {p}Y(x')dx' - / {f} y” dx' ■ a-m(x) , 0<x<c 

2rlb 0 ^ c (3 2) 

■ -y' (x) , c<x<®. 

^ » a^(x-x*) 

where ^ I 

O'^ 

It should be noted that y(x) for 0 < x < c and y(x) for c < x < * are 
Che unknowns. When the jet flap is absent (C^. "0) we obtain Che problem 
solved in section 2.3, whereas when the approaching stream is uniform 
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we obtain the problem solved in ref. 7. The present problem 
appears to be a rather difficult one. The method of solution of the simpler 
problem of ref. 7 with some modifications may still be applicable ‘but this 
needs to be investigated. 


3. RESULTS AND DISCUSSION 


The effect of proxiraity of a surface of discontinuity on the list of an 
airfoil is shown in figure 7. Results from reference 4 are included in this 
figure for comparison. Glauert (ref. 4) represented the airfoil by a single 
vortex distribution. This the only reason for Glauert' s results not 
agreeing with the present ones. For large values of h/c, as can be seen 
from the figure, both methods give the same results. This trend is to be 
expected. At lower values of h/c, the present results are more accurate. 
For values of h/c very close to zero, the results are, however, not 
reliable. 

Figures 8 and 9 illustrate the effect of a uniform jet on the lift of a 
flat plate airfoil. Figure 8 shows that as the ratio of the jet width to 
the airfoil chord (h/c) increases, the lift ratio (L/Lq) increases and 
reaches a value of 1.0 asymptotically. The lift of the airfoil depends on 
the location of the airfoil in the jet. Figure 9 shows that the lift is 
maximum when the airfoil is on the centerline of the jet and decreases as 
the airfoil is moved away from the centerline. This is obviously a result 
of linearization. 

Changes in the lift and pitching moment due to a nonuniform stream and 
wall effects were reported by Ting and Liu (ref. 5). Figure 10 shows the 
results for the example case of reference 5, obtained by the present analy- 
sis and compared with those of reference 5. The two results should have 
been identical. The small differences that can be noticed in this figure 
are due to a small error that had crept in the results of reference S. When 
this was corrected, their results were identical with the present ones. 

The effect of four surfaces of discontinuity on the lift and pitching 
moment of a flate plate airfoil is illustrated in figure 11. The airfoil is 
assumed to be located on the axis of the central stream. The widths of the 
middle three streams are the same. Even though the velocity profile is 
assumed to be symmetric in this example, the computer program developed to 
solve this problem is quite general and accepts different values for the 
velocity in the streams adjacent to the middle stream increases, the lift 
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(and the pitching moment) ratio of the airfoil increases approaching the 
value of 1.0; this increase can be looked upon as the result of an increase 
in effective jet width. 

The effect of a nonuniform stream of a smooth velocity profile is given 
in figures 12, 13, and 14. Figure 12 shows the effect on the lift of a flat 
plate airfoil due to a jet of a Gaussian velocity profile in a uniform 
stream. In figure 12a, a uniform approach stream corresponds to a ■ 0, and 
hence the value L/l^ at this point is unity. For values of 'a* greater 
than zero, the stream has a jet like velocity porfile and the low velocity 
Streams above and below the airfoil cause the lift ratio to drop. For 
values of 'a' less than zero (but greater than -1.0) the stream has a wake 
like profile. In this case the high velocity stream in the neighborhood of 
the airfoil causes an increase in the airfoil lift ratio. 

Figure 12b shows the effect of the spread of the jet on the airfoil 
lift ratio. For small values of the parameter (d/c) the airfoil lift 
ratio is small. As (d/c) increases, the lift ratio increases and reaches 
the value of 1.0 asymptotically as (d/c) tends to infinity. 

Lifting chars'* ter is tics of a Joukowski airfoil in a nonuniform stream 
have been studied by Chow et al. (ref. 1) by solving Euler's equations. 

Some of their results are given in figures 13 and 14 for comparison with the 
results of the present analysis. It may be recalled that the present analy- 
sis is a linear potential flow analysis. 

As the first example following reference 1, the upstream nonuniform 
velocity profile is represented by 

U(y) ■ [ 1+a exp ( ~(y*-y )^/d^)] (33) 

where 'a' is the ratio between maximim excess velocity and the velocity U 
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(at y * °* ) . Parameter d represents the spread of the velocity nonunifono- 
ity and y^ represents the vertical location of the jet centerline with re- 
spect to the airfoil. Figure 13 illustrates the effect of varying the 
location of the airfoil on the lift of a thin Joukowski airfoil (camber ■ 

5X) at zero angle of attack for three different values of the jet spread 
(d/c). In all these cases the present results show that the lift is maximum 


when the airfoil is located on the centerline of the jet. The minerical 
results of Euler's equations included in figure 13 for comparison show that 
the maximum lift occurs when the airfoil is slightly below the jet center- 
line. Generally there is no satisfactory agreement between the two results 
in figure 13. 

As the second example, the lifting characteristics of the Joukowski 
airfoil in a shear layer between two parallel streams is studied. For anal- 
ysis, the upstren velocity profile is chosen to be 

U(y) ■ [ 1+a tanh ((y-y^)/d)] (34) 

The mean velocity Uq is used as the velocity scale. Tbi parameter d rep- 
resents the spread of the shear layer. As before, y denotes the vertical 

s 

location of the airfoil relative to the upstream profile. Figure 14 shows 
the lift coefficient of the Joukowski airfoil placed in such a stream 
obtained by the present analysis together with the numerical results of 
Euler's equation from reference 1. The results show the expected trend, but 
the agreement between the t\ra results is not satisfactory. 

There could be several reasons for the poor agreement between the pre- 
sent results and those of the numerical solution of Euler's equation. The 
present analysis is a linearized potential flow theory. Linearization may 
have introduced some errors. Yet another reason could be due to possible 
errors in the numerical results of reference 1. For example, in reference 
1, a value of 1.315 has been quoted for cc^ for the Joukowski airfoil (c 
■ 1.808) at zero angle of attack in uniform stream. This implies a c^^ of 
0.727. But the linearized theory gives a c^ of only 0.625. Considering 
these differences, the lack of good agreement in results in figures 13 and 
14 is not surprising. 

When a small perturbation approximation is introduced in the vorticity 
transport equation, we arrive at the following linear partial differential 
equation for the perturbation velocity components. 

U(u-v) +v(U) -0 (35) 

y X X y y 

where U ■ U(y) is the undisturbed nonunifora velocity and u(x,y) and 
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v(x,y) are the peiturbation velocity components. Ihe subscripts x and y 
stand for partial differentiation with respect to x and y respectively. 
Since U *■ U(y) is known, the equation is a linear p.d.e. with variable 
coefficients. The linearized boundary condition on the airfoil is as before 
V • Uy', where y' is the slope of the mean camber line. 

The perturbation velocity components u and v s..tisfy the equation 
(U^ - v^) ■ 0. Hence, it is evident that the present potential flow solu- 
tion of linearized Euler's equation (35) only if is small. 


4. CONCLUDING REMARKS 


A linearized potencial flow cheory haa been applied to the problem of a 
thin airfoil in the middle of five atreama with four aurfacea of diaconr 
tinuity. The theory ia then extended to aolve the problem of a thin airfoil 
in a nonuniform atream having a amooth velocity profile. Reaulta have been 
obtained for two examplea - one for a atream having a Gauaaian velocity 
profile and the other having a hyperbolic tangent velocity profile. Theae 
reaulta have u^en compared with more aophiaticated nuaterical reaulta of 
Euler* a equation. Some differences are noticed between theae two reaulta, 
and possible reasons for this poor agreement are indicated. 
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Figure 1. Airfoil Below a Surface of Discontinuity 
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Vortex System tor an Airfoil Below A Surface of 
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Figure 6. Airfoil in a Stream of Smooth Velocity Profile 

















OF POOR QJA. ii , 


d/c=2.0 
d/c=1.0 

U=U«(l+a exp(-(y/d) ^) ) 
c=airfoil chord (=1) 

LQ=Lift in a uniform stream of velocity Uoo(l+a) 
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Figure 12a. Effect of Nonuniform Stream on the Airfoil 

Lift (Variation with the jet excess velocity 
factor a) 
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APPENDIX A 


ORlUir'iAL . 

OF POOR QUALilY 


DERIVATION GF THE RECURRENCE RELATION FOR THE STRENCTH OF IMAGE VDRTICES 


Conaider the re 1 at Iona (16a) to (16d) 


^g+2 “ *^a+2 " ®2 "^1-a 


J ■ r - a r 
"^a+l a+1 1 -a 


Bn C ■ K - o K 
'’0 a N 0 -a-1 


B_iK-L,-a,L 
* a-1 a-1 -1 -a-2 


I • 0, for a > 2 

a 


Jl "0. Tg • r , K-1 ■ 0 , 


L ■ 0, for a < -2 

a 


Ue can rewrite (16b) aa 


r , - B , j , ♦ o ,r 
a+1 1 a+1 1 -a 


The relation (16a) with a replaced by (a-1) givea 


^a+1 “ '^a+1 ‘ “2 "^2-a 


Since 1^^^ ■ 0 for a > 1, thia can be written aa 


“ ■ Vi ‘“2 ■' 2 -. ‘ ' 


or J ,• a, J, 
a+1 2 2-a 


Usiag (A3) in (A2) we obtain 


ORiCil:;^L I 
OF POOH Q,,;.;,;, 


r +<*,r ■ for s > 1 


1 2-8 1 -8 


tihich can be rewritten with s replaced by (1-s) as 


To “S,j, +o,r , 

2-8 1 2-8 1 8-1 


Combining (A5) and (A4) we obtain 


r -a,r ,)+o,r 

8+1 22-8 1 8-1 1 -8 


or r “o,r. ,^+o^«^. a, or 

8 1 -(s-l) 2 -( 8-3) 1 2 8-2 


Next since L ® 0 for s < -2, the relation (16d) reduces to 
$ 

0-1 K«_i = -“-1 L_g_2 for 8 < -1 


Also replacing 8 by -(s+1) in (16d), we obtain 


0-1 *^-(8+2) ’ ^-(s+2) '^-l ^8-1 

Using (A7) in (A8) and noting L. » 0 for s < -1, we get 

(s-17 


K =-o K, -V for 8 < -1 
s-l -1 -(s+2) 


or K ■ - a K , for s < -2, 
s -1 -(s+3) 


Next replacing s by -(s+1) in (16c), we obtain 


^0 "-(s+1) ’ *^-(8+1) ■ “0 ^ 


If this relation is used in (16c), then we find 


^-(s+l) * - Tg + 00 K, 


(A4) 


(A5) 


(A6) 


(A7) 


(A8) 


(A9) 


(AlO) 


(All) 
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ORiGtKi'- < 

OF poor^ qualhY 

which can also be written with s replaced by -(s+3) as follows; 

“0 ^s+2 " ~ ^-(s+3) * *^-(8+3) 

On using (All) and (A12) in (A9), we arrive at 

^8 * “ “-1 ^-(s+3) " “0 ^-(s+1) ■ “0 ^8+2 

(A6) and (A13) are the required recurrence relations for Fs for s > 2 and 
s < - 2, respectively. It can be easily shown with Fq * r that Ti = 

^ -1 “ F. Other values of F'g can be determined using the above 
recurrence relations. 
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